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, Abstract. We perform three-dimensional numerical simulations of surface flows on the companion star in a semi- 

■ detached binary system and construct the corresponding Doppler maps. The mass ratio of the binary system, 

q = M2/M1, considered here is q — 0.05, 0.33, 0.5, 1, 2, and 3. For all cases, we obtain the H-, L1-, and L2-eddies, 
as found in our previous work, and confirm that the flow pattern does not heavily depend on the mass ratio. 
We propose that this kind of problem be dubbed "stellar meteorology." The Doppler maps at the position of the 
companion show a structure tilted towards clockwise direction and presenting deviations from the critical Roche 
surface due to the Ll-eddy and the L2-eddy on the companion star. 

We apply our results to the Galactic supersoft X-ray source RX J0019. 8+2156 and try to attribute the low radial 
velocity component of the emission lines of He II A4686 observed recently to the irradiated spot on the surface of 
the companion rather than that of the white dwarf or the accretion disc. Based on the comparison between the 
observations and our constructed Doppler map, we estimate the mass of the companion star in RX J0019. 8+2156 
to be ~ 2Mq assuming the mass of the white dwarf star to be around O.6A/0. 



Key words, binaries: close - star: evolution - stars: mass-loss - stars: emission-line - stars: individual (RX 
J0019. 8+2156) - stars: winds, outflows - X-rays: stars - accretion, accretion discs 

1. Introduction However, Davey & Smith ( 1992 ) demonstrated the sur- 

face mapping method of the companion surface in cata- 
A semi-detached binary system consists of a Roche-lobe- clysmic vari ables ( CVs) using Na I line absorption. Dhillon 
filling companion star and a (compact) accreting object. & Watson (J2000J) demonstrated a technique of imaging 
Part of the gas of the companion star is gravitationally companion stars of CVs using Roche tomography. Besides, 
attracted by the primary, flows through the LI point to- Doppler map technique also can potentially give some 
wards the primary star and, if the latter is smaller than trace of the surface flow (as considered in this paper), 
the circularisation radius, forms an accretion disc around 

Oka et al. (2002) performed three-dimensional simu- 
lations of surface flows on the companion star and pre- 
Many studies on the accretion disc itself exist, but un- dieted the existence of three kinds of eddy associated with 
til recently, only a very few studies considered the surface a high/low pressure on the companion star: the H-, L1-, 
flow of the companion, except for the pioneering work by and L2-eddies as schematically shown in Fig.l. The no- 
Lubow & Shu l|1975fl . A possible reason for the lack of tations H, LI and L2 denote the high pressure around a 
studies on the companion is that there had been no effec- (north) pole, the low pressure around the LI point and 
tive way to observe the surface of the companion star as the low pressure at the opposite side to the LI point, re- 
well as the flow on it. spectively. 
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Fig. 1. A schematic diagram of the surface flow of a com- 
panion star in a semi-detached binary system: H, LI and 
L2 denote eddies associated with a high/low pressure, re- 
spectively. Typical velocity vectors are shown with num- 
bers. 



In a rotating fluid, the pressure gradient force balances 
the Coriolis force, and the wind blows thus along isobaric 
lines. On the Earth, such a flow is called the geostrophic 
wind, while on a star it was called the astrostrophic wind 
by Lubow & Shu H1975[) . Therefore, a clockwise-rotating 
eddy forms around a high pressure region, while a counter- 
clockwise rotating eddy appears around a low pressure 
region on the northern hemisphere in our counterclockwise 
rotating system (viewed from the north) . We dub this kind 
of problem "stellar meteorology." 

With the observational progress mentioned above as 
well as theoretical ones in hand, we are now in the state 
of investigating the flow on the surface of the companion 
star in some more detail. 

Oka et al. (2002) considered only the mass ratio of 
unity. In reality, however, the range of mass ratio is wide, 
and it is thus interesting to investigate the dependence 
of the flow pattern on the mass ratio, q — M2/M1. In the 
present study, we consider six cases for the mass ratio: q — 
0.05,0.33,0.5,1,2, and 3. As described above, a Doppler 
map can potentially give some traces of the surface flow, 
and so we construct Doppler maps in the present study as 
well. 

1.1. Observational implication of the surface flow on 
the companion star 

Davey & Smith (1992) found an asymmetric irradiation 
pattern in the secondary star of dwarf novae (e.g. IP Peg); 
the trailing hemisphere has stronger Na I line absorption 
than the leading hemisphere (see also Dhillon & Watson 
2000). They considered that this asymmetry could be due 
to the circulation current on the companion star (but see 
Watson et al.EDDSJ- 



As was stated above, Oka et al. 820 02|) found three 
types of eddies on the companion star. They considered 
that the H-eddy as well as the Ll-eddy could transport 
the heat due to the irradiation to the leading hemisphere, 
and argued that these currents naturally explained the 
asymmetry of the irradiation pattern. Note, however, that 
magnetic CVs, on the other hand, show the irradiation to 
be on the trailing hemisphere (Davey & Smith 1992 1996 
Watson et al. I2003J) . To understand the irradiation pat- 
tern of magnetic CVs, we need another mechanism, which 
could lead to the effective heating of the trailing hemi- 
sphere. 

1.2. Supersoft X-ray sources 

Supersoft X-ray sources (SSS) consist, like CVs, of a Roche 
lobe filling companion star and a white dwarf. X-rays are 
believed to be produced by a steady nuclear burning oc- 
curring on the surface of the white dwarf. This requires a 
high mass accretion rate (<^lO -7 M0/yr) that is only pos- 
sible with a high mass companion star (van den Heuvel et 

ai.rr™. 

In our galaxy four SSSs are known (see e.g. Gansicke 
et aL EfflHty : one of which is RX J0019.8+2156. This object 
and its 12.2 mag optical counterpart were discovered by 
Beuermann et al. I|1995[) . Recent observations revealed a 
low radial velocity component to the He II A4686 emission 
line (e.g. Becker et al. H33SI Cowley et al. lTMSl Deufel et 
al. Cnnni McGrath et al.EEDJ. These authors considered 
that the He II A4686 emission line was emitted from the 
white dwarf or from the inner part of the accretion disc. 
With this assumption they estimated the mass of the com- 
panion to be smaller than the white dwarf. However, this 
estimate certainly contradicts the massive companion the- 
ory developed by van den Heuvel et al. (1992). 

Hachisu & Kato (2003), on the other hand, considered 
that the low radial velocity of the He II A4686 emission line 
was emitted from the inner hemisphere of the companion 
star. The inner hemisphere is irradiated by the compan- 
ion and becomes hot enough to emit the He II A4686 line. 
The Doppler map of He II A4686 observed by Deufel et 
al. ( 1999) showed that the emission line peak is located 
around V x ^ 100 km/s and Vy~ 100 km/s (Note that the 
plus sign of the above velocity is due to our choice of the 
coordinate frame) . This emission line peak does not corre- 
spond to the region of the white dwarf or the inner part of 
the accretion disc. With the irradiated companion model 
Hachisu & Kato (2003) argued that such low radial veloc- 
ity emission lines were naturally explained by the massive 
companion star in a system with a mass ratio of q^3 — 4 
and a white dwarf of 0.6M Q . 

Matsumoto & Mennickent (2000) constructed the 
Doppler map of He II A4686 of RX J0925. 5-4758 which 
is another known Galactic SSS and also found a peak lo- 
cated around V* ~ 100 km/s and V y ~ 100 km/s (in our 
velocity sign): i.e., the low radial velocity region. The low 
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Table 1. Mass ratio, computational region, and number 
of grid points of models 



Mass ratio Computed region Number of grid points 

O05 1.0 x f.O x 0.5 101 x 101 x 51 

0.33 1.0 x 1.0 x 0.5 101 x 101 x 51 

0.5 1.1 x 1.0 x 0.5 111 x 101 x 51 

1 1.2 x 1.0 x 0.5 121 x 101 x 51 

2 1.4 x 1.2 x 0.6 141 x 121 x 61 

3 1.4 x 1.2 x 0.6 141 x 121 x 61 



velocity of the He II A4686 emission line seems therefore 
to be a common feature of supersoft X-ray sources. 

Hachisu & Kato (2003) tried to interpret the low ra- 
dial velocity component of the He II A4686 emission line 
in the Doppler map of RX J0019. 8+2156. However, their 
numerical model was two-dimensional, and could there- 
fore not handle flows on the surface of the companion. A 
three-dimensional calculation is necessary in this respect, 
and this is the aim of the present research. 

This paper is organized as follows. In Sect. 2 we de- 
scribe the assumptions of the model and the numerical 
method. In Sect. 3 we show the flow patterns and con- 
struct Doppler maps. In Sect. 4 we discuss the problem of 
the supersoft X-ray sources. A discussion and a summary 
are given in Sect. 5. 

2. Assumptions and Numerical Method 

2.1. Assumptions 

The binary system considered here consists of a primary 
compact star with mass M\ and a companion star with 
mass M<z. The system rotates counter-clockwise viewed 
from the north. The mass ratio of the binary system is 
defined by q = M 2 /Mi. 

We solve the three-dimensional Euler equations in the 
rotating frame of the binary. The equation of state consid- 
ered here is that of an ideal gas characterized by a specific 
heat ratio 7. Complex effects such as physical/turbulent 
viscosity (except numerical one), magnetic fields, and ir- 
radiation from the accretion disc or the primary star are 
neglected. 

The equations are normalized with length by the sep- 
aration A between the centres of two stars, and time by 
l/f2, where fi denotes the orbital angular velocity of the 
binary system and, therefore, the orbital period equals 27r. 
The density on the inner (numerical) boundary is taken to 
be unity. The gravitational constant G is eliminated using 
the above normalization. 

2.2. Numerical method 

We use a Cartesian coordinate system. The origin of the 
coordinate is located at the centre of the primary star. 
The x-axis coincides with the line joining the centres of 
the two stars. We define the x — y plane as the orbital 
plane, so that the z-axis is perpendicular to the orbital 



plane and is oriented in the same direction as the angular 
momentum vector of the orbital rotation. Note that we put 
the mass-donor star in the negative x domain following 
our convention; other authors seem to prefer to put it in 
the positive x domain. Therefore, our sign of velocity is 
different from the convention of many other authors. 

The computational region is a rectangular box, and 
its size and the number of grid points depend on the cases 
as described below. We assume symmetry about the or- 
bital plane, and thus only the upper half of the region is 
computed. We use the simplified flux splitting (SFS) tech- 
nique proposed by Jyounouchi et al. (QJ593} and Shima & 
Jyounouchi (1994) as a Riemann solver, and the MUSCL- 
type technique is used as an interpolation. The method is 
the same as in Makita et al. (2000), Matsuda et al. {2000), 
Fujiwara et al. I|2()()lj> and Oka et al. {5002). The spatial 
and temporal accuracy are kept at second-order levels. 

2.3. Boundary and initial conditions 

We apply the same boundary conditions as in Oka et al. 
(2002). The inner boundary is assumed to be an equipo- 
tential surface slightly smaller than the critical Roche lobe 
(~ 90% of the critical Roche lobe radius) of the compan- 
ion. Note that this inner boundary does not necessarily 
mean the real surface of the companion, but is chosen for 
numerical reasons. 

The inside of the companion star is filled with a gas 
with zero velocity, density po = 1, and sound speed Co, 
where cq is a free parameter. We mainly examine the case 
of Co = 0.05. The gas can naturally flow out from the inner 
boundary if the gas pressure below the inner boundary is 
higher than above it. 

The outside of the outer boundary is assumed to be 
always filled by a gas of fixed velocity, 0, density, pi, and 
sound speed c\, where p\ and ci are parameters that are 
taken to be 10 -5 and vTO, respectively. We note that for 
the calculation of the accretion discs with q=l, 2, and 3, 
the density of the initial ambient gas (pi — —5) is some- 
what large, and therefore affects the accretion discs. In 
order to reduce the effect of the initial ambient gas we use 
a much lower ambient gas density of p\ = — 7 when q is 1, 

2, and 3. 

At the initial time t = 0, the entire region, except the 
inside of the inner boundary, is occupied by a gas with 
velocity 0, p\, and c\. We follow the time evolution until 
a steady state is reached. We describe only steady states 
in the following. 

3. Numerical Results 

3.1. Surface flow 

In this section we describe the pattern of the surface flow 
of a companion star. The cases of q = 0.05,0.33,0.5,1,2, 
and 3 are investigated. The coordinate sizes and number 
of grid points for each case are listed in Table 1. Note that 
the case of q = 1 was already studied by Oka et al. ( 2002). 
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Fig. 2. Equidensity surface (logp = —2) of the compan- 
ion star and the streamlines starting from the equidensity 
surface of logp = —2.5; A mass ratio of q = 0.33 and 
a specific heat ratio of 7 = 5/3, are adopted. The top 
panel is viewed from the north, while the bottom panel 
is viewed from the opposite side of the Ll-point. H-, L1-, 
and L2-eddies can be seen. 

For the calculation of the surface flow we adopt 7 = 
5/3, i.e. we assume adiabatic expansion of the gas in the 
companion's atmosphere. For the calculation of the accre- 
tion discs, on the other hand, we use 7 = 1.01 (almost 
isothermal) to mimic radiation cooling, necessary to form 
an accretion disc. 

Figure 2 shows an equidensity surface of the compan- 
ion star and streamlines originating from another equiden- 
sity surface, for q = 0.33 and cq = 0.05. We can see the 
H-, L1-, and L2-eddies as described above; we refer to Oka 
et al. I|2002|l for the mechanism of generation of these ed- 
dies. It is remarkable that the flow pattern is essentially 
the same as the case of q = 1 in Oka et al. ( 2002). These 
eddies are the manifestation of the astrostrophic wind pre- 
dicted by Lubow & Shu 

Figure 3 depicts the equidensity surface of the accre- 
tion disc with c = 0.05. We can see a pair of spiral shocks 




Fig. 3. Equidensity surface (logp = —3.1) of the accretion 
disc is shown; q = 0.33 and 7 = 1.01 are adopted in this 
figure. The size of the whole computed region is 2.0 x 
1.0 x 0.5 (including companion star), and the number of 
grid points is 201 x 101 x 51. Note, however, that the region 
to the right of the LI point is shown to exhibit only the 
accretion disc. We may observe a pair of spiral shocks and 
a shock associated with the LI stream. 

as well as a shock associated with the LI stream. Figures 
4 and 5 represent the case of q = 2 with Co = 0.05. 
Compared with Figs. 2 and 3, we observe that the flow 
pattern is essentially the same as for the lower q case. 

Figure 6 shows the result of the extreme mass ratio, 
q = 0.05 with Co = 0.02. In this case the gravitational 
potential due to the companion is so shallow that a small 
fraction of the gas escapes through the L2 point. We do not 
discuss here this case in detail. We can conclude that the 
flow pattern does not heavily depend on the mass ratio, 
q, at least for q greater or of the order of 1. 

We note an interesting possibility. As can be seen in 
Figs. 3 and 5, the equidensity surface rises abruptly be- 
hind spiral shocks and is thus making cliffs. These cliffs 
may be irradiated by the radiation emitted from the cen- 
tral star/the central region of the accretion disc, and thus 
produce spiral patterns in the Doppler map. We stress 
here that the spiral shock is a natural consequence of su- 
personic flow in a non-axisymmetric gravitational poten- 
tial. 

3.2. Construction of Doppler maps 

To construct observable Doppler maps, information on the 
ionisation state of atoms in the photosphere is required, 
and the ionisation states are mainly determined by the 
temperature. The temperature in the photosphere of the 
companion and the accretion disc is affected by the ir- 
radiation from the primary and the central part of the 
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Fig. 4. Same as in Fig. 2 but for a mass ratio of q = 2. 




Fig. 5. Same as in Fig. 3 but here the equidensity surface 
of log p = —4.5 is shown for a mass ratio of q = 2. 



Fig. 6. Surface flow on the companion star in a binary 
with q = 0.05; equidensity surface (logp = —2.5) of 
the companion star and the streamlines starting from the 
equidensity surface of logp = — 3. Note that a small frac- 
tion of the gas flows through the L2 point as well as the 
LI point because of the shallow gravitational potential. 



accretion disc. Since, in the present study, we do not take 
the irradiation effect into account, we cannot construct 
real observable Doppler maps but, instead, we use the fol- 
lowing convention. 

Instead of identifying a photosphere, we use an 
equidensity surface as a possible candidate for the photo- 
sphere. The horizontal component of gas velocity, V* and 
V y , of each computational cell at that surface is mapped 
onto a,Vx — V y plane: the Doppler map. We tested various 
values of the density and found that the general trend was 
not affected by a particular choice of the density. 

Our construction method of the Doppler map does not 
reflect the mechanism of line emission, and therefore, only 
a part of our constructed Doppler map can be observed. 
For instance the outer (non-irradiated) side of the Roche 
surface may not be observable, if appropriate lines are not 
emitted from there. 

Another problem in constructing a Doppler map is the 
choice of 7. We use the case of 7 = 5/3 to calculate the 
contribution from the companion, while we use the case 
of 7 = 1.01 for the accretion disc. Thus, we combine two 
different simulations into one Doppler map. Therefore, we 
use two different values of density on an equidensity sur- 
face: a surface of the companion star and one for the ac- 
cretion disc. Since our main aim is to investigate Doppler 
maps caused by the surface flow on the companion, the 
contribution from the accretion disc should here only be 
considered as a reference. 
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Fig. 7. Doppler map constructed from the horizontal com- 
ponents of velocity, V x and V y , at the cquidcnsity surface 
of logp = —2.5 for the companion star (denoted by the 
cross) and \ogp = —3.9 for the accretion disc (denoted 
by the plus). In this figure q = 1. Note that the follow- 
ing figures are placed upside down in order to compare our 
results with those of other authors. The critical Roche sur- 
face and the ballistic orbit of the LI stream are also shown 
for the purpose of reference. The upper and the lower dot 
show the centre of mass of the companion and that of the 
primary, while the middle dot (at the LI point) denotes 
the common centre of gravity. 



Figure 7 shows the constructed Doppler map of the 
cquidcnsity surface of log p = —2.5 for the companion star 
and logp = —3.9 for the accretion disc. The mass ratio is 
q = 1 in this case. Note that all the Doppler maps shown 
in the figures are placed upside-down to accommodate the 
more common convention used by observers. The critical 
Roche lobe of the companion star and the ballistic orbit 
of the LI stream are also shown for reference purposes. 
The crosses show the contribution from the companion, 
while the plus signs show the contribution from the ac- 
cretion disc, which is just for reference as was described 
above and we do not give a detailed discussion about this 
contribution in the present paper. The contribution due to 
the accretion disc collides with that of the surface flow in 
the Doppler map. However, this does not necessary mean 
that the flow in the accretion disc really collides with the 
companion in the configuration space. 

The most important feature in our calculated Doppler 
maps is the asymmetry about the Vy-axis. This is due to 
the surface flow of the companion. If there was no flow 
on the companion surface, all crosses should reside within 
the critical Roche surface. However, this is not the case. 




Fig. 8. The surface of the companion is divided into four 
quadrants as shown in the right bottom sub-panel. We can 
examine the contribution from each surface region. The 
number in each Doppler map corresponds to the quadrant 
in the sub-panel. Note that the same data is used as in 
Fig. 7. 



We can observe that the crosses deviate to the upper right 
and to the lower left. 

In order to see this more clearly, Fig. 8 shows the con- 
tribution in the Doppler map from each quadrant of the 
surface region, described by the numbers in Fig. 8. The 
number in each Doppler map corresponds to the quadrant 
labeled by the same number in the right bottom sub-panel. 
One may safely speculate that the asymmetries seen in 
quadrant 1 and in quadrant 4 are due to the flow in the 
Ll-eddy. The deviations seen in quadrants 2 and 3 are due 
to the L2-eddy, although this portion of the Doppler map 
is unobservable because of lack of irradiation. 

We may summarize the characteristics of the Doppler 
map due to the companion as follows: 

1. The velocity vector denoted by 1 in Fig. 1 has a con- 
tribution of positive V x , V y . This produces the asym- 
metries. 

2. The velocity vector denoted by 2 has a contribution of 
positive V x and negative V y . 

3. The velocity vector denoted by 6 has a contribution of 
negative V x and V y . 

4. The velocity vector denoted by 7 has a contribution of 
positive V x and negative V y . 

5. Velocity vectors associated with the H-eddy, i.e. 3, 4 
and 5, have a tendency to shrink the cross dots in the 
companion towards the centre of the companion. 



Oka et al.: Numerical Simulation of the Surface Flow on a Companion Star II 



7 



4. Supersoft X-ray sources 

We apply our results to the Galactic supersoft X-ray 
source: RX J0019. 8+2156 and try to attribute the low 
radial velocity component of the emission lines of He II 
A4686 in the Doppler map observed by Deufel et al. I|1999f) 
to the irradiated spot on the surface of the companion. 

Figure 9 shows our constructed Doppler map with q = 
3, as a model for the RX J0019. 8+2156. We assume the 
system velocity to be AQ OI b — 330km/s. The equidensity 
surface is chosen with log p — —2.5 for the companion and 
logp = —4.4 for the accretion disc. The region where a 
strong He II A4686 emission line peak around V x ^ 100 
km/s and 100 km/s (with our velocity convention) is 
observed by Deufel et al. (11999J1 . is denoted by a dotted 
circle. 

From Fig. 9, we may conclude that the case of q ~ 3 
can naturally explain the low radial velocity component 
of the emission lines of He II A4686. This is the same 
result as obtained by Hachisu & Kato (2003). The mass 
ratio, q ~ 3, means that the mass of the companion star 
is larger than that of the white dwarf. Considering that 
a mass above ~ 0.8M Q for the WD is ruled out by its 
relatively low X-ray flux of ~ 0.4 x 10 37 ergs/s (Beuermann 
et al. 119951 Meyer-Hofmeister et al. H998fl . we may assume 
the mass of the white dwarf to be ~ 0.6M Q . Therefore, 
the mass of the companion star can be estimated to be 
- 2M Q . 

Finally, as was discussed earlier, the emission peak as 
denoted by the dotted circle is due to the Ll-eddy, and if 
so, we may say that the Doppler map obtained by Deufel 
et al. lfHT99) exhibits the trace of the surface flow on the 
companion star. 

5. Summary 

We investigate the surface flow of the companion star 
in a semi-detached binary system and construct model 
Doppler maps. Six cases with mass ratios of q = 0.05, 0.33, 
0.5, 1, 2, and 3 are considered. We apply our results to the 
Galactic supersoft X-ray source RX J0019.8+2156. The 
results are summarized as follows: 
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1. 



We obtain the H-, L1-, and L2-eddies, which were 
found by Oka et al. I|2002(l . for all mass ratios. We 
confirm that the flow pattern does not heavily depend 
on the mass ratio. 

We construct model Doppler maps, and we find an 
asymmetric velocity distribution, due to the surface 
flow. 

The companion mass of RX J0019. 8+2156 is estimated 
to be around 2Mq. This massive companion is consis- 
tent with the theoretical model proposed by van den 
Heuvel et al. 



4. We argue that the Doppler map of supersoft X-ray 
sources shows indications of the surface flow on the 
companion star. 
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